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ABSTRACT 

We have discovered a remarkable warm (130 — 160 K) molecular hydrogen tail with a H2 mass 
of approximately 4 x 10 7 M extending 20 kpc from a cluster spiral galaxy, ESO 137-001, in Abell 
3627. At least half of this gas is lost permanently to the intracluster medium, as the tail extends 
beyond the tidal radius of the galaxy. We also detect a hot (400 — 550 K) component in the tail that is 
approximately 1% of the mass. The large H2 line to IR continuum luminosity ratio in the tail indicates 
that star formation is not a major excitation source and that the gas is possibly shock-heated. This 
discovery confirms that the galaxy is currently undergoing ram-pressure stripping, as also indicated 
by its X-ray and Ha tails found previously. We estimate the galaxy is losing its warm H2 gas at a 
rate of ~ 2 — 3 M Q yr _1 . The true mass loss rate is likely higher if we account for cold molecular 
gas and atomic gas. We predict that the galaxy will lose most of its gas in a single pass through the 
core and place a strong upper limit on the ram-pressure timescale of 1 Gyr. We also study the star- 
forming properties of the galaxy and its tail. We identify most of the previously discovered external 
Ha sources within the tail in our 8/im data but not in our 3.6/im data; IRS spectroscopy of the region 
containing these Ha sources also reveals aromatic features typically associated with star formation. 
From the positions of these HII regions, it appears that star formation is not occurring throughout the 
molecular hydrogen tail but only immediately downstream of the galaxy. Some of these HII regions 
lie outside the tidal radius of the galaxy, indicating that ram-pressure stripping can be a source of 
intracluster stars. 

Subject headings: infrared: galaxies — galaxies: individual (ESO 137-001) — galaxies: clusters: indi- 
vidual (Abell 3627) — galaxies: evolution — galaxies: ISM 



1. INTRODUCTION 

The morphologies, gas content, and star-forming char- 
acteristics of galaxies vary drastically from the cluster to 
the field environment (see Boselli & Gavazzi (2006) and 
the references therein) . There are a range of possibilities 
regarding the observed gas deficiency and reduced star 
formation of cluster spirals. One possible explanation is 
that gas-rich field galaxies are transformed as they fall 
into a cluster, during which process they are stripped of 
their gas and their star formation is halted. A variety of 
processes, such as ram-pressure stripping (Gunn & Gott 
1972) and tidal interactions between galaxies (Moore et 
al. 1996) and with the cluster potential (Henriksen & 
Byrd 1996), could be at play. However, a complete 
picture of the transformation of galaxies as they enter 
the cluster environment is still missing because we do 
not understand the degree to which cluster specific pro- 
cesses affect the evolution of the galaxies within them. 
Several studies suggest that the relative importance of 
ram-pressure stripping and tidal interactions depends 
on the mass of the galaxy in question (Cortese et al. 
2007; Haines et al. 2006, 2008). Theoretical work shows 
that the transformation probably starts in group envi- 
ronments (Kawata & Mulchaey 2008), but gas loss also 
occurs from gas-only stripping (ram-pressure) events in 
the cluster environment itself (Tonnesen et al. 2007). The 
relative roles of these environments are ambiguous and 
would be clarified by an improved understanding of the 
effects of ram-pressure stripping. Truncated Ha disks 
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and HI tails on galaxies with undisturbed stellar disks 
provide evidence for this process (Koopmann & Kenney 
2004; Chung et al. 2007). However, our understanding 
of the physical nature of the stripping from such systems 
will be improved by additional direct measures of the 
interaction with the intracluster medium. 

Spitzer can offer a unique view by observing the inter- 
action between the ICM and galactic interstellar medium 
(ISM) directly. Appleton et al. (2006) serendipitously 
discovered powerful molecular hydrogen emission ema- 
nating from the intragroup medium (IGM) of Stefan's 
Quintet. This emission is remarkable because the H2 
dominates strongly over other line emission (Cluver et 
al. 2010) and the IR continuum is weak. The emission 
is thought to arise from a high velocity (~ 1000 km s _1 ) 
collision between an intruding galaxy and a tidal arm 
(Guillard et al. 2009). Appleton et al. (2006) place con- 
straints on the temperature components of the excited 
gas and its column density by measuring the pure rota- 
tional H 2 transitions S(0) thru S(5). The observed H 2 
lines are shock-excited because the spectrum is domi- 
nated purely by molecular hydrogen lines and the sur- 
face luminosity of the H2 is ten times that from the 
X-ray emission associated with the shock- front. They 
show that there are at least two temperature compo- 
nents in the molecular gas: a warm and dense component 
(~ 200K) and a hot and sparse component (~ 700K). A 
more recent work by Cluver et al. (2010) that uses Spitzer 
data covers a significant portion of the Stefan's Quintet 
and confirms and expands on the conclusions of Appleton 
et al. (2006). 

Observing these molecular hydrogen lines (if they ex- 
ist) around galaxies falling into clusters would probe the 



2 



gas-stripping processes and would provide quantitative 
estimates of the amount of gas being stripped. An excel- 
lent candidate for this test is ESO 137-001, a 0.2L* SBc 
galaxy in the Norma cluster (Abell 3627). ESO 137-001 
is a remarkable galaxy with coaligned X-ray and Ha tails 
that extend approximately 70 and 40 kpes, respectively 
(Sun ct al. 2006, 2007, 2010). Its projected distance from 
the center of the Norma cluster is only 0.28 Mpc. Its low 
radial velocity (77 km s -1 ) with respect to the cluster as 
a whole means most of its motion is in the plane of the 
sky. Moreover, deep Ha imaging reveals bright extra- 
planar knots along the Ha tail that are most likely very 
large star- forming regions (Sun et al. 2007), suggesting 
the stripped gas may be forming stars. Deep Rp-band 
imaging also shows low surface brightness streams that 
are aligned with the X-ray tail but have a larger width 
(Woudt et al. 2008). All evidence suggests this galaxy is 
currently undergoing ram-pressure stripping. 

In this paper we report the discovery with Spitzer of a 
warm molecular hydrogen tail that trails ESO 137-001. 
Detailed Spitzer observations of this galaxy and its tail 
expand on our understanding of ram pressure stripping 
and permit the estimation of the resulting mass flow. 
The paper is structured as follows. In §2, we discuss the 
details of our observations and data reduction. In §3, 
we present the results of our analyses. In §4 we discuss 
the implications of our results. Finally in §5, we list our 
conclusions. For computing distances, we adopt the con- 
cordance cosmological model = 0.73, O m = 0.27, 
and .ffo = 71 km s" 1 Mpc" 1 ). We adopt redshifts pro- 
vided by the NASA Extragalactic Database of 0.0157 
and 0.01544 for the cluster and galaxy, respectively. This 
yields a luminosity distance of 67.1 Mpc to the cluster 
and an angular scale of 1" = 0.316 kpc. All reported 
errors in this work are quoted at the la level. 

2. OBSERVATIONS AND DATA REDUCTION 

In this paper we detect and characterize warm H2 emis- 
sion and resolve the morphology of the stripped gas in 
the tail of ESO 137-001, and study the star- forming prop- 
erties of the galaxy and the extraplanar HII regions. To 
this end, we carried out IRAC imaging and IRS spec- 
tral mapping observations. Our IRAC imaging allows 
us to resolve star-forming regions and study the stellar 
properties of the galaxy. IRS spectral maps provide de- 
tailed information (at a coarser spatial resolution) about 
the physical properties of the ISM such as the strength 
of the aromatic features, molecular hydrogen lines, and 
fine-structure lines. The data were taken as part of the 
Spitzer GTO program 50213 (PI: G. Rieke), as discussed 
below. 

2.1. IRAC Imaging 

We observed ESO 137-001 on May 10, 2008 in all four 
(3.6, 4.5, 5.7, and 8/zm) IRAC channels (Fazio et al. 
2004). Frame times were kept short to prevent the sat- 
uration of the array from Milky Way foreground stars 
because the galaxy is located at a low galactic lati- 
tude (b <~ —7°). We obtained a total on source in- 
tegration time of 260 s in each of the channels. We 
constructed IRAC mosaics of the galaxy with MOPEX 
(version 18.1.5) using the Basic Calibrated Data (BCD) 
from the Spitzer pipeline (version S17.2) as input and 



following the Spitzer Science Center (SSC) IRAC reduc- 
tion cookbook. Bad pixel rejection was achieved through 
dithering and median filtering. 

2.2. IRS Spectral Mapping 

The observations of ESO 137-001 used both the SL 
and LL IRS modules that together span 5.3 to 38.0 fim 
(Houck et al. 2004). The full wavelength range allows us 
to observe several ground vibrational state H 2 rotational 
lines, specifically the z^=0-0 S(0) thru S(7) transitions. 
We used the most recent pipeline reduced data (version 
S18.7) and CUBISM (version 1.7) spectral map cube 
construction software (Smith et al. 2007) for all of our 
analyses to take advantage of the most up-to-data data 
calibrations. Our first set of IRS observations on May 
2, 2008 was designed so each sky position was visited 
by two different IRS pixels — the minimum required to 
produce a good spectral map. This was achieved by step- 
ping the SL and LL slits in the perpendicular direction 
by 1.85" and 5.1", respectively, such that the spectral 
map was centered about the galaxy. Moreover, during 
the time of our observations the LL slits were aligned 
to include the galaxy tail. For our second set of IRS 
SL-only observations on April 6, 2009, we used the same 
mapping technique, but covered a region immediately 
west of the galaxy (presumably on the leeward side of 
the ICM wind). We used the outrigger off-source point- 
ings from both LL and SL observations to measure the 
background. Figure 1 shows the final configuration of 
the IRS slit positions for both sets of observations. Fi- 
nal on-source integration times per pixel for our first set 
of IRS observations were 235 s, 235 s, 2265 s, and 880 
s for the SL1, SL2, LL1, and LL2 spectroscopic cubes, 
respectively. The final per pixel integration time for our 
second set of IRS observations was 940 s for each of the 
SL1 and SL2 spectroscopic cubes. 

Background subtraction was carried out using CU- 
BISM, which averages the frames of all off-source 
pointings associated with a given spectrograph using 
a min/max clip, and subtracts it from the on-source 
frames. If there is a bright source within one of the off- 
source pointings, the averaging technique of all off-source 
pointings diminishes the contribution of the source. The 
final background frame was visually inspected; the LL1, 
SL1, and SL2 (for both datasets) background frames 
showed no obvious signs of bright source contamination. 
However, the LL2 background revealed a few faint con- 
tinuum sources that were on the south east side of the 
galaxy. Fortunately, these sources do not contaminate 
the regions from which we extract our spectra in this pa- 
per, and we do not attempt to mitigate this issue. All SL 
and LL spectrographs suffer from significant bad pixel 
issues, which we addressed through a manual process. 
We visually inspected the spectral cube for any artifacts 
along the direction of the slit steps. We manually masked 
the offending pixel that produced the artifact (identifying 
it using the backtracking feature in CUBISM), rebuilt 
the spectral cube, and verified the artifact was removed. 

The spectra from the different spectrograph orders 
were stitched together and placed on a single wavelength 
scale without any need for rescaling the flux of any indi- 
vidual order. Our processing of the data produces fluxes 
in the orders that are consistent with each other in the 
regions of overlap. As a check, we compared the IRS 
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fluxes with an independent 8/im IRAC flux measure- 
ment. For the spectral map centered around the galaxy, 
we found the 8/im flux (derived from a weighted aver- 
age using the IRAC 8/zm filter transmission as weights) 
of the final reduced spectra of a 22" diameter aperture 
centered about the galaxy agreed with the IRAC 8^m 
extended-source corrected flux to within 15%. For the 
first IRS dataset, the sizes of the final spectral maps are 
67" x 30", 183" x 41", and 20" x 28" for SL-only, LL- 
only, and both SL and LL maps, respectively. For the 
second IRS dataset, the sizes of the final spectral maps 
were 59" x 16.7", 39" x 18.4" (parallelogram shaped) for 
SL-only, and both SL and LL maps. When the LL data 
were combined with the SL data, we interpolated the LL 
data onto the SL data's pixel grid. CUBISM generated 
error cubes associated with each spectral map. We used 
these error values in our spectral extractions unless oth- 
erwise noted. They do not include systematic errors, for 
example due to calibration errors or fringing. 

3. RESULTS 

3.1. Conditions in the Warm Molecular Hydrogen Tail 

In the lower right panel of Figure 2, we show the IRS 
image taken at the rest-frame 17.035^™ H 2 — S(l) 
transition. We generated this image by summing up the 
flux from 17.2 — 17 Afim, and show only the emission that 
is 3cr above the background. We find a remarkable tail 
emanating from ESO137-001 that extends all the way to 
the edge of the spectral map. This tail, aligned with 
both the Ha and X-ray tails discussed by Sun et al. 
(2006, 2007), is at least 20 kpc in length and approxi- 
mately 3 — 5 kpc in width. The virtually undisturbed 
stellar morphology of ESO 137-001 (Woudt et al. 2008) 
plus the unusual state of the gas in the tail eliminate 
the possibility that it is a tidal feature (see Sections 4.1 
and 4.2). Instead, our observations support previous sug- 
gestions that the tail arises from ram-pressure stripping 
(Sun et al. 2006, 2007). We present spectra of this tail in 
Figures 3 and 4, which were generated by summing up 
the spectra within the LL-only (shown by the large or- 
ange and smaller red rectangles in Figure 2) and SL/LL 
(shown by the red polygon in Figure 2) regions, respec- 
tively. The large LL-only extraction aperture is called 
the LL-only full tail region, while the smaller LL-only 
aperture is called the LL-only far tail region. The solid 
angles (f2) of the SL/LL, LL-only full tail, and LL-only 
far tail are 9.07 x 10" 9 , 2.38 x 10 Js , and 1.09 x 10~ 8 sr" 1 , 
respectively. The LL-only far tail region was chosen to 
be uncontaminatcd by star formation (see Section 3.2.1) 
and outside the tidal radius of the galaxy. The spectra 
show strong detections of molecular hydrogen rotational 
(H2 v = — 0) lines in these regions: the 5.511 fim 
S(7), 9.665 \im S(3), 12.279 fim S(2), 17.035/xm H 2 0-0 
S(l), and 28.22fim H 2 0-0 S(0) transitions. Several fine 
structure lines and aromatic features were also observed 
within these apertures. 

3.1.1. Molecular Hydrogen Tail Properties 

We calculate the physical properties of the stripped 
molecular hydrogen by assuming an optically-thin, single 
(or double) temperature gas model. We first compute the 
column density of each line, N s ^, using the following 
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where flux is determined using PAHFIT (version 1.2; 
Smith et al. 2007), which fits for both the continuum 
and lines, Ag^ is the Einstein A-coefficient for a given 
transition obtained from Turner et al. (1977), vs(i) 1S 
the frequency of the transition (values are tabulated in 
Black & van Dishoeck (1987)), and f2 is the solid angle of 
the aperture used. The errors for the line fluxes are ob- 
tained from la errors reported by the fitting algorithm. 
We do not include any extinction correction in the fit of 
the lines because it is likely to be low in the tail region 
outside of the galaxy. We present the excitation diagram 
{\n.{N S u\/ g u ) versus T u ) from the observed ground-state 
rotational transitions in the LL-only regions and SL/LL 
region in Figure 5. g u and T u represent the degeneracy 
and transition energy expressed in terms of the tempera- 
ture of the upper level of the transition. The degeneracy 
g for molecular hydrogen can be expressed as (2 J + 1) 
for even/para states (i.e. J = even) and 3(2 J + 1) for 
odd/ortho states (i.e. J = odd). We have firm detec- 
tions of the S(0) and S(l) lines in the LL-only region 
and S(0), S(l), S(2), S(3), and S(7) lines in the SL/LL 
region. From the excitation diagram, we can infer that 
there are multiple temperature components in the molec- 
ular hydrogen tail within the SL/LL region. 

Due to the lack of electrical dipole moment in the hy- 
drogen molecules, the only allowed radiative transitions 
between rotational states require AJ = ±2. Therefore, 
ortho and para states are radiatively decoupled. The 
ortho/para ratio (OPR) of molecular hydrogen in lo- 
cal thermodynamic equilibrium (LTE) varies from ~ 1.5 
at 100 K to - 3 above 200^ (Burton et al. 1992). 
The lack of radiative coupling requires collisions between 
molecules (and/or a catalyst) to reach LTE. The warm 
component of the molecular hydrogen has temperatures 
that typically fall within this range of temperatures. 
Consequently, the OPR may not conform with the LTE 
value if the hydrogen molecules formed at temperatures 
higher than 200 K and have now cooled to temperatures 
lower than 200 K without sufficient density to reach LTE. 
The OPR is then frozen to the value when the gas was 
last in LTE. We carry out a test outlined by Roussel 
et al. (2007) to check if the observed population is not 
in LTE. They characterize any deviation from LTE with 
the O-P-Rjngh T parameter, which is the OPR at high gas 
temperatures. At LTE, O-P-R^ighT is by definition 3. 
Roussel et al. (2007) state that if the following condition 
is satisfied the gas is consistent with a population in LTE: 
T(5(0) - 5(1)) < T(S(0) - S*(2)) < T(5(l) - S*(2)) < 
T(5(l) - 5(3)) < T(S(2) - 5(3)) where T(S(i) - S(j)) 
is the temperature derived from S(i) and S(J) line ra- 
tios. The line ratios of lower excitation energy lines probe 
components at lower temperatures in a gas with a dis- 
tribution of temperatures. Temperatures derived from 
T(5(0) - 5(2)) and T(5(l) - 5(3)) are independent of 
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whereas the rest depend on it. The temper- 



atures derived from our lines satisfy the above equality 
and are consistent with being in LTE. This does not mean 
the population is definitively in LTE, as it is possible to 
have slightly lower values of the OPi? m gj 1T ratio and 
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still satisfy the condition. However, we assume LTE in 
our fits of the molecular hydrogen population. 

To quantify the physical properties of the gas, we first 
determined the temperature and total column density 
(N tot ) of the LL-only regions by fitting the S(0) and S(l) 
lines using Equations 2 and 3: 



ln(N(T u )/g u ) = In(JV ) - 



T 



(2) 



where N(T U ) is the column density of a transition, g u 
is the degeneracy of the transition, Nq is the column 
density of the J = state, T u is the energy (in Kelvin) 
of the upper energy level of the transition, and T is the 
temperature of the gas. 
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where N tot is the total column density of the gas, j is the 
rotational state, Tj is the energy of that state. T = 
K, and T\ = 1 70.6 K for molecular hydrogen. Equation 
3 was summed up to the H2 0-0 J — 25 state, which cor- 
responds to an excitation energy of 37728 K, well above 
the observed temperatures. Given the limited informa- 
tion, we fitted a single temperature of 157 ±3 K and a 
total column density of 2.46+°;23 x 1Ql9 cm ~ 2 in the LL- 
only full tail region. The errors in the temperature and 
column density were determined through Monte Carlo 
runs. 10,000 fitting trials were calculated where the trial 
line fluxes were drawn from a gaussian distribution with 
a mean and a full-width-half-max (FWHM) that were 
defined by the measured line flux and line flux error, re- 
spectively. We derive a total H2 mass within the tail 
of 4.1lt°;^ x 10 7 Mq. In Figure 5, the blue line is the 
best fit single-temperature model. For the LL-only far 
tail region, we obtain a temperature of 149 ± 6K and 
a total column density of 2.90t°;43 x 1Ql9 cm~ 2 . This 
corresponds to a total H 2 mass of 2.15l ; 32 x 10 7 M©. 

Next, we characterized the molecular hydrogen pop- 
ulation in the SL/LL region. As shown by Figure 5, 
the additional information requires a model with multi- 
ple components with different temperatures. We chose 
the simplest model consisting of two temperature com- 
ponents as given in Equation 4. 



\n{N{T u )/g u ) = ln(iV 0;1 )+ 
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where A^i is the column density of the J = state 
of the first component, T\ and T 2 are the temperatures 
of the first and second components, respectively, and 
/2,i = N n _ 2 /N ,i where iV ,2 is the column density of 
the J = state of the second component. We only fit 
the S(0), S(l), S(2), and S(3) lines. We did not include 
the S(7) line in our fits for two reasons: 1.) we would 
require a third temperature component to account ade- 
quately for this line, but we do not have sufficient line 
detections to constrain a third component; and 2.) this 
line is most affected by UV fluorescence, which is a re- 
sult of star formation. Our spectral maps showed some 
evidence for a patchy distribution of the S(7) line emis- 
sion, though no clear correlation was observed with star- 
forming regions. The best fit model (shown in black in 
Figure 5) provides temperature estimates of 1251^ K 



and 472^5q K for the two components. The total column 
density for each component was calculated using Equa- 
tion 3, producing a value of 5.9t 4 ' 2 x 10 19 cm~ 2 for the 

warm and 5.9+^6 x 10 17 cm ~ 2 for the hot components. 
These values correspond to masses of 3.6±jj;? x 10 7 M© 
and 3.6^ 2 'g x 10 5 M© for the warm and hot components, 
respectively. The results of all fits are tabulated in Table 
3. These results agree well with the hydrodynamical sim- 
ulations of ram stripping by Roediger & Briiggen (2008). 
They predict that massive stripped tails can be produced 
even 0.5 to 1 Mpc from the cluster core (the projected 
distance of ESO 137-001 is 0.3 Mpc) with column densi- 
ties of ~ 10 19 cm~ 2 . 

We also calculate the total H2 line luminosity in each 
extraction aperture. To make this measurement, we must 
make a small extrapolation to include the luminosity 
from ground state rotational lines that are at shorter 
wavelengths than those already detected. We carry this 
out using the following equation: 



L(H 2 
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where is the solid angle of the aperture, Da is the angu- 
lar size distance to the galaxy, N s ^ is the column density 
of the molecules at the upper energy state for the S(i) 
transition, and A and v are the Einstein A-coefficient and 
frequency associated with that transition. The column 
densities are calculated from values obtained from the ex- 
citation diagram fit. For the SL/LL region, we obtain a 
luminosity of 1.59x 10 40 ergs s _1 , while for the LL full tail 
and far tail regions, we obtain 1.36 x 10 40 and 0.59 x 10 40 
ergs s _1 , respectively. Even though the LL full tail re- 
gion covers the entire tail (i.e. larger solid angle than 
the SL/LL region), its H2 luminosity is lower than the 
SL/LL region because of the significant emission com- 
ing from the hot component within the SL/LL region. If 
there is a hot component associated with either of the 
LL regions, then the stated luminosities are lower limits. 
Furthermore, there may be a third even hotter compo- 
nent inside the SL/LL region as revealed by the poor fit 
of the S(7) line, which suggests that even the luminosity 
derived for the SL/LL region is a lower limit. When com- 
pared to other wavelengths, we obtain luminosity ratios 
of L(H 2 )/L(Ha) ~ 1 and L(H 2 )/L (X-ray) - 0.1 where 
the L(Ha) and L(X-ray) are obtained from Sun et al. 
(2007) and Sun et al. (2006), respectively. It is clear H 2 
is a significant coolant. These ratios are only approxi- 
mate, as the apertures over which the comparisons are 
made are not completely matched; the H2 aperture is 
smaller than those used in the Ha and X-ray measure- 
ments, so its cooling effect may be underestimated. 

3.1.2. Comparison of Molecular Hydrogen Gas Mass 

We compare our results with the other mass estimates 
of the tail derived from observations in the Ha-band (Sun 
et al. 2006), and the X-ray (Sun et al. 2007). To make 
these comparisons, we compute the average gas surface 
density to mitigate issues arising from mismatched aper- 
tures. The gas surface density derived from our molecu- 
lar hydrogen study ranges from 0.4 — 0.9 M© pc~ 2 . Our 
measurement of molecular hydrogen mass is relatively 
robust because it is reasonably independent of the tail 
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geometry and does not require a free-parameter /, the 
filling factor. Nevertheless, there are systematic issues 
in our mass estimate: 1 .) the IRS slit coverage does not 
cover the full axial extent of the tail, which causes us to 
miss some of the H 2 emission at the edges; and 2.) it is 
unclear what fraction of the molecular gas the warm H2 
emission represents. Roussel et al. (2007) in their sur- 
vey of warm molecular hydrogen in SINGS galaxies find 
that for galaxies with CO detections (43 galaxies in to- 
tal) warm molecular hydrogen (T > 100 K), on average, 
represents only ~ 10% of the cold molecular hydrogen. 
There is a very broad range of warm to cold gas ratios 
in the SINGS sample ranging from 0.6 — 45% with a 
scatter for the full sample of a <~ 10%. It is unclear if 
we can apply the fraction found in SINGS galaxies to 
explain the molecular hydrogen that is being stripped 
when the underlying formation and heating mechanism 
may be different. Preliminary CO measurements of the 
shocked region of Stefan's Quintet suggest a higher frac- 
tion (~ 50%) of warm to cold molecular gas (Guillard 
2010). However, Stefan's Quintet is an extreme envi- 
ronment with exceptionally strong molecular hydrogen 
emission, so the observed fraction may not be represen- 
tative of conditions in the warm H2 tail of ESO 137-001. 

Both the Ha and X-ray estimates rely on additional 
assumptions that we do not require. Both assume a 
cylindrical distribution for the gas to estimate an av- 
erage electron density, which is then used to arrive at 
the mass of the radiating gas. Moreover, both estimates 
depend strongly on /, which is unknown. Taking these 
effects into consideration, we compute a surface density 
of hydrogen by dividing the mass estimate by the angular 
extent of the tail and obtain ~ 2.5 /s M Q pc~ 2 for both 
the Ha and X-ray tails. For the X-ray mass value, we 
use the tail mass from Sun et al. (2006). This is approx- 
imately a factor of 3 — 6 higher than our estimate if we 
assume / = 1, but agrees well with ours for the plausible 
value / = 0.03 — 0.1, assuming all three tails have similar 
masses. 

The exact physical mechanism that excites the 
stripped gas to produce emission at these different wave- 
lengths is not entirely clear, making it difficult to predict 
the different states of the gas (see Section 4.1 for a dis- 
cussion). One needs to carry out a full census of the 
stripped gas through CO and HI observations to obtain 
a complete picture. Cold atomic gas measurements of 
the galaxy have been carried out, but were only able to 
place an upper limit on the HI mass of 10 9 M Q (Sun et 
al. 2007). If we assume the typical values found in the 
SINGS sample, cold H2 may be a significant component 
of the stripped gas, and its mass is critical for calculat- 
ing an accurate mass loss rate for the galaxy (see Section 
4.2.3). 

3.1.3. Fine Structure Line Emission in Tail 

We also fit several common fine structure lines us- 
ing PAHFIT in the LL-only and SL/LL-only regions to 
study possible excitation mechanisms for the observed 
molecular hydrogen tail. The lines are highlighted in Fig- 
ures 3 and 4, and the line fluxes are tabulated in Table 2 
for all three spectral extraction regions. We detect [Nell] 
(12.8/xm), [Nelll] (15.6/xm), [SIII] (18.7 and 33.5/xm), 
and [Sill] (34.8/xm). We also observe an emission line in 
the SL/LL region at rest 26 fim, which is likely a blend 



of the [OIV] 25.89/jto and the [Fell] 25.99/jm lines. It 
is possible to distcntangle which of the two lines domi- 
nates the blend. The [OIV] line and the [Fell] line have 
vastly different excitation potentials, 56 and 7.9 eV, re- 
spectively. It is likely that this line detection is domi- 
nated by the [Fell] line as the low [NeIII]/[NeII] ratio of 
0.2 suggests the gas is not highly excited. Cluver et al. 
(2010) draw similar conclusions where they too observe 
a 26/jm line in Stefan's Quintet. Further progress can be 
made by comparing the [SIII] (33.5/xm)/[SiII] ratio. Dale 
et al. (2006) carried out an IRS survey of a variety of dif- 
ferent types of active galactic nuclei (AGN) and nuclear 
and extranuclear HII regions. They found a strong cor- 
relation between the type of object and the [SIII] /[Sill] 
flux ratios. HII dominated spectra typically have higher 
[SIII] /[Sill] ratios than AGNs and therefore this ratio 
can be a good discriminator of the type of energy source 
exciting the gas. We evaluate where each of the three 
spectral extraction regions fall in this correlation. The 
SL/LL region has a [SIII] /[Sill] ratio of 0.75, which is in 
a region that is largely dominated by HII regions. The 
same is true for the LL-only full tail, which has a ratio of 
0.67. However, the LL-only far tail region yields a much 
lower ratio of 0.32, which is in an area that is dominated 
by AGNs. This is another indication that the gas in this 
region may be excited by a mechanism other than star 
formation, as discussed in Section 4.1. 

3.2. Star Formation in ESO 137-001 and its Tail 

The Ha image clearly shows that the galaxy is actively 
forming stars in the nuclear region and to a lesser extent 
along the tail. Our infrared data offer us a complemen- 
tary look at the star formation activity. We detect 8/j,m 
excess sources, which are clearly associated with the Ha 
sources found by Sun et al. (2007). We measure the 
SFR of the galaxy, focusing on the nuclear region where 
we see the majority of the star formation as indicated 
by Ha imagery and we have full wavelength coverage 
across 5 — 36/im. Finally, we address two questions: 1.) 
Is ESO 137-001 peculiar in any way in its star-forming 
properties? 2.) How does star formation progress in the 
molecular hydrogen tail? 

3.2.1. 8 ixm Dust Emission 

We present the results of our IRAC imaging in Figure 
2. The Ha image of ESO137-001 was kindly provided 
by M. Sun and has been discussed in detail by Sun et 
al. (2007). There are striking differences between the 
3.6 \im and 8 /im images. The 3.6/im data shows what 
appears to be a normal spiral galaxy with significant con- 
tamination from foreground Milky Way stars. The star- 
subtracted, contrast-enhanced, deep Rp-band image of 
Woudt et al. (2008), reproduced in Figure 1, shows this 
galaxy more clearly and reveals only minor low surface 
brightness perturbations that extend in the southwest 
direction. However, the 8/zm image shows a completely 
different view. We see a cometary structure centered 
about the galaxy within which are several point-sources. 
Excess 8 fim flux often arises from warm dust rich in 
aromatic hydrocarbons, associated with star formation 
(Calzetti et al. 2007). 

To understand the nature of the 8/zm sources, we first 
remove the 8)im stellar flux to determine the emission 
from the aromatics. We achieve this by first convolving 
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the 3.6/im data to degrade the image quality to match 
the 8/im resolution using kernels provided by Gordon et 
al. (2008). Then we subtract the convolved 3.6 /im image 
after registering it to the 8/im image and multiplying 
it by an arbitrary scaling (~ 0.3) to ensure the bright 
stars are properly subtracted. The final subtracted im- 
age is presented in Figure 6. The kernel is imperfect 
because the wings of point-sources are over-subtracted; 
however, we are able to discern 8/im excess clearly. We 
determine the position of the 8/zm excess sources using 
SExtractor (Bertin & Arnouts 1996). We first identify 
sources that have a minimum of 4 pixels detected more 
than 1.5ct above the local background. We apply an ad- 
ditional 5c cut to these sources using the source flux and 
flux error generated by SExtractor to generate the source 
catalogue. We visually inspect the final source catalogue 
and remove the sources that are associated with detector 
and stellar continuum subtraction artifacts. This method 
identifies all of the sources that can be identified by eye. 
We present the locations of the sources within 2' of the 
galaxy in Table 4 along with their distances from the cen- 
ter of the galaxy (a =16:13:27.30 and S = -60:45:50.5) 
obtained from the NASA Extragalactic Database. 

A comparison with the Ha emission (see Figure 6) 
shows that the structures and point-sources seen at 8/zm 
are also present in the Ha image. The 8 /tm excess 
source positions are overlaid on the Ha image to indi- 
cate the correspondence. There are a few bright pixels 
surrounded by dark annuli seen at 8/zm, which are in fact 
residuals from the stellar continuum subtraction. When 
we do a careful comparison of the two images almost ev- 
ery source at 8/im within 1' (approximately within the 
tidal radius) of the galaxy has a clear Ha counterpart 
with one sole exception, a point source located southeast 
of the galaxy, also located on the side of the galaxy op- 
posing the tail. It is unclear what this source is, and it 
is not very likely to be associated with the galaxy. The 
few fainter Ha sources that are not seen in 8 /zm can 
attributed to the reduced sensitivity in the 8 /zm band 
due to its lower angular resolution. Overall, the source 
matching strongly suggests that these 8/im excess sources 
are associated with the galaxy and are relatively young 
HII regions. This confirms Sun et al. (2007) 's hypothesis 
that these are star-forming regions (which they arrived 
at from their analyses of the B-band, I-band, and Ha 
data), and possible sites for intracluster star formation. 

3.2.2. Measurement of Star Formation Rate in the 

Galaxy 

We calculate the star formation rate (SFR) in the cen- 
tral parts of the galaxy by first estimating the MIPS 
24/zm (Rieke et al. 2004) flux density from our IRS spec- 
trum multiplied by the MIPS relative response curve. 
We find a MIPS 24/zm flux density of 45 ± 2 mJy in 
the nuclear extraction region, which contains almost all 
of the 24 /zm flux (see Figure 8). The uncertainty in 
the flux was determined by computing the standard de- 
viation of the residual of a continuum fit from 20 to 30 
/im, excluding the S(0) line. This is necessary because 
the CUBISM software underestimates the errors in that 
wavelength range and small ripples associated with fring- 
ing are present in the continuum. Hereafter, all MIPS 
24/im flux density measurements are carried out in this 
manner. 



Using Rieke et al. (2009) 's SFR-to-24/im calibration 
produces a SFR value for the galaxy of 0.6 M Q yr" 1 . 
We note, however, our galaxy's 24/im luminosity falls a 
factor of 2 lower than the valid range of calibration of 
the relation in Rieke et al. (2009). By using their cali- 
bration, we likely underestimate the total (UV+IR) SFR 
somewhat because the UV leakage from the star-forming 
regions will be higher for lower 24 /zm luminosity galax- 
ies. According to the Figure 7 of Buat et al. (2007), the 
UV leakage can be about 15% for our galaxy. After tak- 
ing this into account, we deduce the total SFR (UV+IR) 
of ESO 137-001 to be 0.7 M Q yr" 1 . This value is domi- 
nated by the systematic uncertainty of the SFR-to-24/zm 
calibration ~ 0.2 dex (Rieke et al. 2009). 

Our SFR value is a factor of 2 lower than that mea- 
sured by Sun et al. (2007). Sun et al. (2007) used an 
extinction-corrected Ha measurement to estimate a SFR 
of 2 M yr -1 , which translates to 1.3 M Q yr -1 for the 
Kroupa IMF used by Rieke et al. (2009). There is a 
mismatch in the aperture used by Sun et al. (2007) and 
our work for the measurement of the SFR. We correct 
for this using the Ha data provided by M. Sun to mea- 
sure the relative difference of total Ha flux between the 
two apertures. The difference is very small; the Ha SFR 
measurement is reduced by 5% when the 24/im aperture 
is used yielding a Ha SFR of 1.2 M Q yr -1 . No statistical 
uncertainty for this measurement was quoted by Sun et 
al. (2007); however, this measurement carries a system- 
atic error of <~ 0.15 dex (Kennicutt 1998). Sun et al. 
(2007) applied a 1.4 mag extinction correction estimated 
from the H7/H/3 line ratio. Reconciling the Ha with the 
IR-dctcrmincd SFR, however, would require that there 
is 0.8 mag of extinction at Ha. To probe this inconsis- 
tency, we use the Kroupa IMF SFR calibration given by 
Kennicutt et al. (2009) that uses both Ha (not extinc- 
tion corrected) and 24/im luminosity. Sun et al. (2007) 
provide an extinction corrected L(Ha) within the region 
they calculate the SFR. After removing the extinction 
correction and applying the aperture correction, we ob- 
tain a L(Ha) of 6.5x 10 40 ergs s _1 . Using this value along 
with the 24/zm luminosity yields a SFR of 0.7 M Q yr -1 , 
which is consistent with our original 24/im SFR measure- 
ment. This suggests that the extinction correction used 
by Sun et al. (2007) was probably overestimated. 

We now compare the SFR with the stellar mass of ESO 
137-001 to see if it has an abnormal SFR. We first esti- 
mate the galaxy's stellar mass from its K-band magni- 
tude. We obtained a K-band apparent magnitude of 11.6 
from a catalogue of galaxies in Abell 3627 (Skclton et al. 
2009), which corresponds to a Mk = —22.5. This yields a 
stellar mass of - 1.4 x 10 10 M (Bell & de Jong 2000). A 
typical SFR for a cluster member of this mass is <~ 5 M Q 
yr -1 , with a board distribution extending down to the 
rate in ESO 137-001 (Vulcani et al. 2010). Thus, there 
arc no obvious anomalies in the star-forming properties 
of the galaxy, although it is less active than is typical. 

To explore the galaxy's star formation properties fur- 
ther, we look to its IR spectrum (see Figure 7). The 
galaxy unsurprisingly shows strong aromatic features and 
an infrared continuum, both clear indicators of on-going 
star formation. We compare the spectrum with tem- 
plates provided by Smith et al. (2007) and find it is well- 
matched with a star-forming galaxy template as shown 
by Figure 7. We also checked to see if there are any pecu- 
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liarities in the aromatic features of ESO 137-001. We fit 
the nuclear spectrum of the galaxy using PAHFIT (ver- 
sion 1.2; Smith et al. 2007) with extinction correction, 
and the best fit is shown in Figure 7. We tabulate the 
strength of various aromatic features in Table 5. The de- 
rived flux ratios L(6.2/im)/L(7.7/im Complex) of 0.3 and 
L(11.3/im Complex) /L(17/um Complex) of 2.1 are consis- 
tent with star-forming galaxies with Hll-dominated nu- 
clei (Smith et al. 2007). 

3.2.3. Star Formation in the Molecular Hydrogen Tail 

The Ha image (Sun et al. 2010) demonstrates that star 
formation is occurring outside the galaxy near the head 
of the tail. To explore the star formation in this region, 
we compare the IRAC 8/im data and the IRS 24/im map 
with the molecular hydrogen S(l) line map. We first 
convolve the 8/im continuum-subtracted image with a 
kernel we created using the 17/im IRS PSF generated 
by stinytim and the 8/xm PSF using the same technique 
outlined above to degrade the Ha image. In Figure 8, we 
show that the S(l) emission is much more extended than 
the 8/xm emission, indicating that star formation only 
proceeds in the part of the tail immediately downstream 
of the galaxy. This is corroborated by our 24/im map 
shown in Figure 8. We measure a MIPS 24/xm-band 
flux of 3.5 ± 0.5 mjy within the SL/LL region where the 
majority of the extraplanar SF is occurring. The flux in 
this region is less than 10% of the value at the galactic 
nucleus indicating that most of the SF is confined to the 
galaxy. Even though we detect 24/im emission in the 
SL/LL region we do not quote a SFR. It is unclear how 
to translate the measured flux in the star formation rate 
as the ICM environment is hostile to dust, and the usual 
SFR calibrations may not hold. We also confirmed the 
lack of MIPS 24/zm-band flux detection (< 1.10 mjy, 
3(J detection limit) in the LL-only far tail region, where 
there is no associated 8/xm excess emission and half of 
the detected warm H 2 lies. This is also corroborated by 
the low [SIII] / [Sill] ratio in this region. We conclude a 
significant portion of the stripped molecular gas may not 
have high enough density to support star formation. 

We also looked for any peculiarities in the aromatic 
features measured in the SL/LL region. We tabu- 
late the fluxes in Table 5. The derived flux ratio 
f(6.2/im)/f(7.7/im Complex) of 0.11 is peculiar for a star 
formation dominated region (Smith et al. 2007). Its value 
is lower than the lowest value found by Smith et al. 
(2007) of 0.2 for Hll-dominated nuclei, which have a typ- 
ical value of <~ 0.3. However, the flux ratio f(11.3/xm 
Complex) /f(17/xm Complex) of 1.9 is not inconsistent 
with the average value observed of Hll-dominated nuclei. 
One possible explanation for the observed discrepancy is 
that aromatic emission at 6.2/im is somehow suppressed. 

4. DISCUSSION 

4.1. Excitation Mechanisms for Molecular Hydrogen 

There are a number of mechanisms to excite molecular 
hydrogen: 1.) UV fluorescence (Black & van Dishoeck 
1987); 2.) X-ray irradiation arising when X-ray photons 
penetrate the molecular gas clouds deeply (Maloney et 
al. 1996); and 3.) shock-heating from high-velocity colli- 
sions. In the last relatively fast J-shock can dis- 
sociate H 2 at the shock front only to have it form later 



on grain surfaces downstream. The lowest J rotational 
lines are collisionally excited in this process (see Hollcn- 
bach & McKee (1989)). Recently, Guillard et al. (2009) 
have presented a sophisticated model of the H 2 emis- 
sion from Stcphan's Quintet. They are able to produce 
a multiphase medium that emits strongly in rotational 
molecular hydrogen lines, and to a lesser extent in X-ray 
and Ha, from a large galactic scale shock caused by the 
collision of an intruding galaxy with a tidal arm. The 
shock- heated gas eventually cools to form molecular H 2 , 
which is the main coolant in the gas, and is excited by low 
velocity, non-dissociative magnetohydrodynamic (MHD) 
shocks that dissipate most of mechanical energy of the 
collision. Their model adequately explains the two tem- 
perature molecular hydrogen component and the high 
L(H 2 )/L(IR) ratio discovered by Appleton et al. (2006). 
This model qualitatively matches our observations, and 
we conclude that the H 2 tail can be shock-excited. 

To investigate further the H 2 excitation mechanism, 
we compare our results with the warm molecular hydro- 
gen inventory of the SINGS star-forming galaxies (Rous- 
sel et al. 2007). We use the molecular hydrogen flux to 
24/im flux ratio as a proxy for L(_ff 2 )/_L(/i?). The ratios 
of molecular hydrogen line flux (measured for the 0-0 
S(0) thru S(3) lines) to MIPS 24/im flux for SINGS star- 
forming galaxies are shown as diamond points in Figure 
9. The errors in each data point shown are the errors in 
the molecular hydrogen line flux. Because they should 
be much smaller, the 24/im flux errors were not quoted 
and therefore they are not included in the error in the 
ratio. We also calculated the same ratios for the SL/LL 
and the LL far tail regions, which are shown by the tri- 
angular and square points respectively. The S(2) and 
S(3) line fluxes for the LL region were calculated by ex- 
trapolating the temperature fit using Equation 5. The 
LL region derived ratio is a lower limit, as we only have 
an upper limit for the 24/im flux. The ratios calculated 
for the tail regions are clearly significant outliers when 
compared to the SINGS sample even in the case of the 
SL/LL region of ESO 137-001 where there is contamina- 
tion from star formation. In fact, the ratios are compara- 
ble to or greater than the most significant outlier in the 
SINGS sample: NGC 4450. Roussel et al. (2007) state 
that their two most significant outliers, NGC 4550 and 
NGC 4579, have significant excess in H 2 emission that 
cannot be explained by emission from PDRs or super- 
nova remnant shocks. They state that shocks produced 
by cloud collisions in combination with X-ray irradiation 
could explain the energetics of the observed H 2 emission. 

We summarize several pieces of evidence that suggest 
that some process other than star formation must be the 
source of the H 2 excitation: 1.) the most compelling ev- 
idence is that approximately 50% of the warm molecular 
hydrogen gas mass, which resides outside the tidal radius 
of the galaxy and within the hot ICM, has no detected 
star formation activity; 2.) even the region within the 
tail with some star formation activity shows unusually 
high luminosity of H 2 for its IR luminosity in comparison 
with the average value of SINGS galaxies; 3.) the H 2 line 
emission is at least 10% of the observed X-ray luminosity 
and approximately similar to the Ha luminosity of the 
tail, which suggests a highly energetic process heating 
the gas; and 4.) the [SIII] / [Silll] ratio in the tail out- 
side the tidal radius is consistent with gas excited by an 
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AGN, possibly suggesting that the energy source may be 
hard. We conclude that the H2 tail is likely stripped gas 
from the galaxy that is either heated by X-rays and/or 
through shocks produced by the galaxy's interaction with 
the ICM. We may have directly detected the molecular 
gas from on-going ram pressure stripping of an infalling 
cluster spiral. This possibility is discussed further in Sec- 
tion 4.2.3. 

4.2. Physical Mechanism Responsible for the Tail 

ESO 137-001 has a relatively normal spiral-like stel- 
lar distribution and asymmetric gas distribution as ev- 
idenced by the X-ray and Ha tails. The existence of 
the molecular hydrogen outside the plane of the galaxy 
and its anomalous excitation mechanism raise several in- 
teresting questions: How did this molecular gas find its 
way out of the galaxy? What does this say about galaxy 
transformation in the cluster environment? How long 
does this phenomenon persist? To answer these ques- 
tions we look at the possible processes that strip galaxies 
of their gas in clusters. We discuss why tidal stripping 
processes are not significant, and why ICM gas-stripping 
is the most likely scenario that explains the observed 
phenomenon. 

4.2.1. Galaxy-Galaxy Interactions 

Galaxy-galaxy interactions can be strong transforma- 
tive processes in the cluster environment. These inter- 
actions can be classified into two broad categories: low 
speed and high speed. Low speed interactions frequently 
lead to significant morphological changes and mergers, 
whereas high speed encounters lead to tidal stripping of 
gas and stars from the outer edges of disks and induce 
nuclear starbursts. 

If virialized, the high velocity dispersion of Abcll 3627 
(a = 925km s _1 ; Woudt et al. 2008) makes it an unlikely 
site for a low velocity interaction between two galaxies. 
However, the case of Abell 3627 is not simple. X-ray 
imaging of the cluster reveals a subgroup near its center, 
which is identified in X-ray residuals after subtracting 
a radially symmetric component. The cluster's diffuse 
ICM emission is elongated with a position angle of ap- 
proximately 130° (Bohringcr et al. 1996). The central 
cD galaxy that lies within the subgroup has a pecu- 
liar velocity of 561 km s _1 with respect to the cluster 
mean velocity indicating that the cluster is undergoing a 
merger at its center (Woudt et al. 2008). However, ESO 
137-001 lies outside of the subgroup. There is also no 
compelling evidence that ESO 137-001 is part of a dy- 
namically cooler subgroup (Woudt et al. 2008) making 
it unlikely for it to have experienced a low velocity en- 
counter with a nearby galaxy. Furthermore, the closest 
galaxy located within 10' with the smallest difference in 
radial velocity with ESO 137-001 is ESO 136-024, which 
has a Av r — 120 km s _1 . One would expect tidally dis- 
rupted material to follow the trajectory of the perturbcr, 
but the location of ESO 136-024 is not aligned with the 
tail. In addition, there is no evidence in the stellar mor- 
phology (Woudt et al. 2008, see Figure 1) for distortions 
in the direction of ESO 136-024. 

Galaxy harassment in the form of multiple fast inter- 
actions in clusters can strip gas and stars from small 
spiral galaxies and turn them into dwarf elliptical or 



dwarf spheroidal galaxies (Moore et al. 1996). Simula- 
tions show that for high surface brightness spiral galaxies 
with stellar disk scale lengths of 3 kpc, a value compa- 
rable to that measured by Sun et al. (2007) for ESO 
137-001, the effect of high speed encounters with massive 
galaxies is likely to be minor (Moore et al. 1999). In their 
simulations of harassing tidal interactions involving high 
surface brightness model galaxies similar to ESO 137- 
001, Moore ct al. (1999) showed that the disk lost only 
1—2% of its stars and largely retained its disk shape while 
its disk scale height increased by a factor of 2 — 4 and 
its disk scale length decreased somewhat. The current 
disk morphology of ESO 137-001 suggests that if it had 
an interaction with a massive galaxy, it is likely to have 
had a minor impact on its stellar distribution. There is 
no physical reason why the gas would be preferentially 
stripped over stars in a tidal interaction. Moreover, other 
simulations show that in tidal interactions, gas typically 
tracks stars in tidal tails whereas ram-pressure can pro- 
duce a significant difference between the stellar and gas 
distribution (Vollmer 2003; Mayer et al. 2006). In fact, 
simulations show that tidal interactions often lead to gas 
inflow and bar formation, and ram-pressure is necessary 
to strip the gas (Mayer et al. 2006). Additionally, there 
is only a single tail, whereas both leading and trailing 
tails are usually present in a tidal interaction (Moore et 
al. 1999). However, simulations show that it may be pos- 
sible for ram pressure to erase tidal features produced by 
a single close encounter and produce long gaseous tails 
trailing the galaxy (Kapferer et al. 2008). Therefore, the 
clear difference between the distribution of stars within 
ESO 137-001 and the gas that trails behind it suggests 
that tidal interaction between galaxies is not the main 
contributing factor for the existence of the tail. 

4.2.2. Galaxy- Cluster Tidal Interactions 

Inside the cluster core, the deep cluster gravitational 
potential can also foster galaxy transformation (Henrik- 
sen & Byrd 1996). When the cluster tidal force exceeds 
the galaxy's own restoring force, gas and stars can be 
tidally stripped from the galaxy and trigger a nuclear 
starburst (Henriksen & Byrd 1996; Cortese et al. 2007). 
We calculate the effect of the cluster potential for our 
galaxy by using the formalism outlined by Henriksen & 
Byrd (1996). The cluster tidal effect can be expressed as 
a radial and a transverse acceleration, a r and a t , respec- 
tively, which are given below: 

a t (r) = GM cl (r) — - 2 (6) 



[R 2 + {r + R) 



and 



a r (r) = GM cl {r) 



(r + Rf 



(7) 



where M c [ is the mass of the cluster, r is the distance of 
the galaxy from the cluster center, and R is the radius 
of the galaxy. We can calculate the significance of the 
cluster tidal force by comparing it to the centripetal ac- 
celeration of the galaxy, a ga i, at the outermost parts of 
its disk: 

_ GM gal 
a ga i - R2 (0) 

where M ga i is the dynamical mass of the galaxy, and R 
is the physical extent of the galaxy. On one hand, if the 
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radial tidal acceleration exceeds the centripetal acceler- 
ation of the galaxy, we would expect the gas and stars 
to be stripped by the cluster's gravitation force. On the 
other hand, if the transverse tidal acceleration exceeds 
the centripetal acceleration, the gas clouds within the 
galaxy will be compressed leading to a burst of star for- 
mation. 

To determine the strength of the cluster tidal acceler- 
ation, we first estimate the cluster dynamical mass pro- 
file. We tested two approaches to determine the mass: 
one that uses the X-ray derived dynamical mass profile 
(Bohringer et al. 1996), and the other that extrapolates 
to smaller radii using an NFW profile (Navarro et al. 
1996) and the dynamical mass estimate determined by 
Woudt et al. (2008) at large radii. To find the best- 
fit NFW profile, we use the dynamical mass values for 
the cluster at 3 different radii (0.67, 1.35, and 2.02 
Mpc) from Woudt et al. (2008). We use only the virial 
theorem (Mvt) and robust virial Mrvt dynamical mass 
estimates at these radii for our fit. Woudt et al. (2008) 
suggest that their third method for mass estimation, pro- 
jected mass estimator, is affected most by substructure, 
which is present in Abell 3627, and produces a mass esti- 
mate that is 50% larger than Mvt- Mvt and Mrvt are 
least affected by substructure and agree to within 5%, 
which is why we fit to the average of these two values. 
Nevertheless, these masses are also likely to be overes- 
timates due to the on-going merger in the cluster. The 
NFW parameters we obtained were r s = 346 kpc and 
c = 6.2. We do not quote errors on these parameters 
because error estimates were not given for the mass esti- 
mates. 

We now compute the relative strengths of the cluster's 
radial acceleration and the galaxy's centripetal accelera- 
tion. The stellar disk is detected out to 40" in radius by 
Sun et al. (2007); this corresponds to a physical distance 
of 12.6 kpc, which we pick to be the radius of the disk. 
We estimate the galaxy's dynamical mass by using the 
K-band velocity-luminosity relation given in Courteau 
et al. (2007). For its K-band luminosity, ESO137-001 
should have a rotational velocity of 110 km s _1 , which 
coupled with its stellar disk size, corresponds to a dynam- 
ical mass of 3.5 x 10 10 M Q . In Figure 10, we show the ratio 
of the cluster's radial and galaxy's centripetal accelera- 
tion for the NFW and X-ray gas mass profile. Even at 
cluster-centric distance of 280 kpc, the minimum likely 
separation of the galaxy from the cluster center, we find 
that the cluster tidal field does not exceed the galaxy's 
centripetal acceleration. We also calculate the galaxy's 
tidal radius to be 15 kpc at this cluster-centric distance 
using the NFW profile for the cluster mass distribution. 
This is consistent with the value derived by Sun et al. 
(2007). Therefore, cluster tidal forces are insufficient to 
cause any appreciable stripping of the stellar and gaseous 
components of the galaxy. 

4.2.3. Hydrodynamic Stripping 

We have shown that the molecular tail emanating from 
ESO 137-001 cannot result solely from galaxy interac- 
tions or from cluster tidal effects. We now explore the 
possibility that it arises from ram pressure stripping. The 
fact that the stellar disk appears relatively unperturbed 
signals that the tail must be generated by a mechanism 
that only affects gas. ICM ram-pressure stripping or vis- 



cous stripping are the most likely candidates. 

To evaluate this possibility, the cluster-centric radius 
within which ram-pressure exceeds the galaxy's restor- 
ing force can be estimated using the Fujita & Nagashima 
(1999) reformulation of the Gunn & Gott (1972) crite- 
rion: 

PlCMVg a i > IllGYjstar^HI 

>2.11 x 10- 11 dyne cm" 2 ( " rot X ( — 

V220kms" 1 / V 10 

X \8 x 10 20 m H cm' 2 J 

where pic'M is the density of the ICM, v ga i is the ve- 
locity of the galaxy with respect to the ICM, v rot is 
the rotational velocity of the galaxy, R is the radius 
of the disk, Eg/ is the surface density of the HI disk, 
and rriH is the mass of the hydrogen atom. Bohringer 
et al. (1996)'s /3-model fit (Cavaliere & Fusco-Femiano 
1978) to their ROSAT X-ray data provides an estimate 
of picm- We use the /3- model parameter values of po = 
4.82 x 10~ 27 g cm" 3 , r x = 184.5 kpc, and /3 = 0.555 to 
compute the ICM density as a function of cluster-centric 
radius. We also use the rotational velocity and radial ex- 
tent of the galaxy given in Section 4.2.2. The HI surface 
density and the 3D velocity of this galaxy are unknown. 
Broeils & Rhee (1997) and Karachentsev et al. (1999) 
have shown that the HI surface density of spiral galaxies 
is fairly constant irrespective of the Hubble type or ro- 
tational velocity. Karachentsev et al. (1999) specifically 
looked at a nearby volume-limited sample of galaxies and 
found the HI surface density had a weak dependence on 
a galaxy's rotational velocity at less than the la level. 
The studies give slightly different mean surface densi- 
ties. Broeils & Rhee (1997) compute an average surface 
density of 3.8 ± 1.1 M pc~ 2 for their sample of spirals 
and irregulars. From the HI surface density - rotational 
velocity relation shown in Karachentsev et al. (1999), we 
estimate a surface density of 6.3 M pc -2 . For our calcu- 
lations we were conservative and chose the higher surface 
density, as that increases the minimum pressure required 
for ram-pressure stripping. For the galaxy velocity, we 
picked the values that bound a galaxy's speed in a viri- 
alized cluster: a (925 km s^ 1 ) and V3a (1602 km s" 1 ) 
where a is the velocity dispersion of the cluster. As an 
additional check we compute the Mach number of the 
galaxy using the X-ray data from Sun et al. (2006). As- 
suming the galaxy is experiencing a shock as it travels 
through the ICM, we can predict the Mach number if 
we know the density and temperature of the head of the 
X-ray tail and the ICM. Sun et al. (2006) quote an av- 
erage electron density and temperature for the head of 
the tail to be 1.8 - 3.5 x lO" 2 /" 1 / 2 cm" 3 and 0.64 kcV, 
respectively, where f is the filling factor. The ICM on the 
other hand has a temperature of 6.3 keV and an electron 
density of 1.4 x 10~ 3 cm~ 3 at the location of the galaxy 
(Sun et al. 2006). If we assume a filling factor of 1, we 
obtain a Mach number that ranges from 1.1 — 1.6. Using 
the ICM sound speed of approximately 800 km s _1 that 
we calculate below, we obtain a galaxy velocity through 
the ICM of 880-1280 km s _1 . This matches closely with 
our assumed value, suggesting that the galaxy cannot be 
moving significantly faster than our predictions. 
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For a galaxy traveling at a velocity of 925 km s -1 , 
ram-pressure stripping will become significant at cluster- 
centric radii less than 0.72 Mpc. However, if the galaxy 
is traveling at 1602 km s _1 the radius increases to 1.4 
Mpc. We show these results in Figure 10 for the two 
velocity cases where we plot the effectiveness of ram- 
pressure in terms of the parameter q ram where q ra m — 
PiCMV 2 gal /2irGY, star Y, H i. When q ram > 1, ram-pressure 
stripping becomes effective. Because ESO 137-001 has a 
very low line-of-sight velocity with respect to the cluster 
(~ 77km s _1 ), its true distance from the cluster center is 
not likely to be much larger than its projected distance. 
Boselli & Gavazzi (2006) state that infalling galaxies typ- 
ically follow highly eccentric orbits. We can estimate the 
cluster-centric distance of the galaxy assuming a roughly 
linear orbit with the following equation: 

r= (11) 

cos(sin {vios/vgaij) 

where r is the cluster-centric distance, R is the projected 
distance from the cluster center, vi os the line-of-sight 
velocity of the galaxy, and v ga i the 3D velocity of the 
galaxy. The two galaxy velocity limits mentioned pre- 
viously put the galaxy at a cluster centric distance of 
approximately 280 kpc. The galaxy is well within the 
minimum radius limits where ram-pressure should be a 
significant transformative process. 

We also consider the case where the galaxy is undergo- 
ing viscous stripping (Nulsen 1982). This can be a signif- 
icant stripping process with the same scaling properties 
as ram-pressure stripping if the gas flow is turbulent. 
For turbulent flow, the kinematics of the galaxy through 
a cluster should produce a Reynolds number > 30. The 
Reynolds number, Re, is given in terms of the mean-free- 
path and sound speed in the ICM, A/cm and cjcm'- 



Re = 2X 



( r 9al \ 
\^ICM ) \CICM 



•Jgal 



CICM - 



A/CM = 11 



kT 

TlCM 



fio- 



10 8 K J \picM/m H 



(12) 
(13) 
(14) 



where T is the temperature of the ICM. We assume an 
ICM temperature of 6.3 kcV (Sun et al. 2006). We cal- 
culate Re for the galaxy located at 280 kpc from the 
cluster center for the two velocity limits. We obtain 
30 < Re < 52 for 925 < v gal < 1602 km s _1 . We there- 
fore expect that the gas flow across the wake of the galaxy 
is turbulent, as also shown in hydrodynamic calculations 
(e.g., Roediger & Briiggen 2008). It is plausible that 
eddies in this part of the tail promote the star forma- 
tion occurring there. We can also estimate the mass loss 
rate following the prescription in Nulsen (1982), which 
is M — itr 2 pvgai where r is the radius of the gas disk 
and p is the density of the ICM. We set the radius of 
the gaseous disk to be the half-width of the 17.035 pm. 
emission at the head of the warm H2 tail, which is ap- 
proximately 3 kpc. This is the best guess at the width 
of the current gas disk width. We obtain a mass loss 
rate that ranges from 0.7 < M < 1.2 M yr _1 , which 
is a factor of 3 lower the value that we obtained from 
our warm molecular hydrogen observations (see Section 



4.2.4 for warm H2 mass loss rate). This indicates that 
ram-pressure and not viscous stripping dominates, but 
viscous stripping may also contribute to the mass loss. 

4.2.4. Mass Loss Rate and Stripping Time Scale 

We calculate an instantaneous warm H 2 mass loss rate 
by dividing the measured warm H2 mass (LL full tail re- 
gion) within the 20 kpc long tail (seen in projection in 
the sky) by the time it would take the galaxy to traverse 
that distance. We adopt the galaxy velocity limits of 925 
km s _1 < Vg a i < 1602 km s _1 used in our ram-pressure 
calculation. The galaxy's velocity is taken to lie in the 
plane of the sky due to its low recessional velocity with 
respect to the cluster. It will take the galaxy approxi- 
mately 1 — 2 x 10 7 years to traverse the length of the tail, 
which corresponds to a warm H 2 mass loss rate of 1.9—3.4 
M© yr _1 . However, it is unclear what fraction of the to- 
tal stripped gas the warm H 2 represents; the mass loss 
rate presented is clearly a lower limit as evidenced by the 
Ha and X-ray tails already detected. In order to make 
an estimate of the total gas mass loss rate, we need mea- 
surements of two important parameters: the total H 2 to 
warm H 2 mass ratio (7// 2 ) and the total gas to H 2 mass 
ratio (jgas)- Accounting for these ratios yields a total 
mass loss rate of (1.9 — 3.4)7// 2 7 gas M Q yr _1 . However, 
both of these 7 values are unknown for this galaxy's tail 
because no CO observations or HI measurements of the 
tail exist. X-ray and Ha measurements do provide some 
constraints of the atomic hydrogen content of the tail, 
but the measured mass depends strongly on the filling 
factor, which is unknown. 

To estimate the true mass loss rate of the galaxy, we 
make a few simplifying assumptions to determine 7 val- 
ues by using measurements made in other systems. 7// 2 
can be bound by cold-to-warm H 2 gas mass ratios deter- 
mined by Roussel et al. (2007) for SINGS galaxies and by 
Guillard (2010) for the shock front in Stefan's Quintet. 
In the SINGS sample, Roussel et al. (2007) on average 
measure a jh 2 of ~ 10, while Guillard (2010) measures a 
7zj 2 of ~ 2 at the strongest portion of the galactic shock 
in Stefan's Quintet. This provides a plausible range of 
2 < 7/j 2 < 10. It is likely that the gas in the shock front 
in Stefan's Quintet is in a much more extreme state than 
in our case, yielding a much smaller ju 2 ratio than what 
is likely for the tail of ESO 137-001. This possibility is 
strengthened by the rather large 7// 2 of <~ 10 measured in 
NGC 4450, which is also thought to have a similar exci- 
tation mechanism as our H 2 tail (Roussel et al. 2007). H 2 
is likely to be disassociated during the process of strip- 
ping, in which case the observed molecular hydrogen may 
form inside the tail after the shock-heated gas cools. In 
any case, we do not expect jh 2 to be a strong function 
of time. 

In comparison, it is plausible that "f gas has changed 
as the stripping has progressed. Initially, the large HI 
envelope of the galaxy would have presented a large cross 
section for stripping. There is even indirect evidence 
that the HI loss rate must have been high in the past. 
To illustrate, we estimate the initial gas mass of ESO 
137-001 prior to its entering the cluster. We do so from 
the stellar mass of the galaxy, multiplied by the typical 
value of the gas to stellar mass ratio for a galaxy of the 
given K-band luminosity, M gas = 0.8 M*, where M* is 
the stellar mass of the galaxy (Bell & dc Jong 2000). 
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The result is that the galaxy should have had an initial 
gas mass of ~ 10 10 M , with an rms scatter of 0.2 dex 
(Bell & de Jong 2000). This mass is about an order of 
magnitude higher than the upper limit to the current 
HI content of ESO 137-001 (Sun et al. 2007), indicating 
that most of the gas has already been stripped. This is 
not surprising, since ram-pressure stripping is thought 
to become important at 3 — 5 times the current cluster- 
centric distance of the galaxy. The time during which 
this stripping has occurred is at most 400 — 700 Myr, 
the time to reach its current position at the range of 
estimates for the velocity of the galaxy from the radius 
where the ram-pressure criterion is met. 

We can also estimate an upper limit for stripping the 
remaining gas. We combine the upper limit on the HI 
mass of 10 9 M with an estimate of the initial H 2 mass 
as 0.1 M gas (the expected value for a Sc galaxy with 10 10 
M of gas (Obreschkow & Rawlings 2009)), or 10 9 M , 
to set an upper limit of 2 x 10 9 M for the current gas 
content (HI plus H 2 ). From the loss rate of warm H 2 and 
the lower limit of ~ 2 on jH2 , there is a lower limit to the 
current mass loss of 3.8 — 6.8 M yr _1 , even assuming 
all the lost gas is in molecular form (that is, ignoring 
the additional mass loss for gas in atomic form). We 
then find an upper limit of 300 — 500 Myr to lose the 
remaining gas. This timescale is similar to the time it 
would take the galaxy to travel from its current position 
to the outskirts of the cluster core, which we take to be 
the core radius given by the X-ray gas distribution, for 
the range of galaxy velocities explored. This means the 
galaxy is currently roughly half way through the process 
of being completely stripped of its gas and will lose most 
of its gas before it finishes its first pass toward and then 
away from the cluster core. 

This also yields an upper limit on the timescale for 
ram-pressure stripping of <~ 1 Gyr, which is the sum of 
the time the galaxy has been experiencing ram-pressure 
stripping and the time it will take to lose all of its 
gas from now, for the range of galaxy velocities ex- 
plored. There are several reasons why the true stripping 
timescale is much shorter: 1.) We assume no molecular 
gas has been stripped in our calculation of the time it 
would take for the galaxy to exhaust its gas from the 
present moment. 2.) The current measured HI mass 
used in our calculation is only an upper limit. 3.) We 
assume no atomic gas is currently being stripped, which 
does not explain the existence of both a X-ray and a Ha 
tail. 4.) We do not account for the impact of star forma- 
tion, which could be a significant consumer of gas. 5.) 
Even though the ram-pressure criterion may be met at 
large cluster-centric distances, the ICM wind may not 
be strong enough to strip the gas permanently from the 
galaxy's halo during the early stages of stripping. This 
is consistent with simulations of ram-pressure stripping, 
which indicate that significant gas stripping happens in 
a relatively quick timescale of a few hundred Myr for a 
galaxy of similar gas mass but 10 times more massive 
(Roediger & Hensler 2005; Kapferer et al. 2009). This 
may suggest that the warm H 2 may only represent a 
small fraction of the total gas stripped. 

4.3. Implications 

The existence of H 2 in the ram-pressure stripped tail 
is puzzling because this is at odds with conventional wis- 



dom that molecular gas is much more tightly bound than 
HI gas, making it far more difficult to strip than the less 
dense atomic gas. This is supported by the observation 
that cluster galaxies with strong HI deficiencies also have 
fairly normal amount of molecular gas content (Bosclli 
& Gavazzi 2006). To determine how molecular hydrogen 
can exist outside the galaxy within the hot intracluster 
medium we carry out a simple calculation of the effect 
of ram-pressure on a typical molecular cloud. At the 
current location of ESO 137-001 where ram-pressure is 
significant, a spherical molecular cloud 10 pc in radius 
and 10 4 M in mass will only be moved by ram-pressure 
a few parsecs in a 1 Gyr. This means that ram-pressure 
cannot move whole clouds out of the galaxy. It is possible 
that molecular clouds may be ablated by the ICM wind 
and stripped gradually. However, the more likely expla- 
nation is that the high velocity dissociative shocks aris- 
ing from the interaction between the ICM and the ISM 
of the galaxy dissociate the molecular hydrogen within 
the galaxy only to have it form on grain surfaces when 
the shocked gas cools (Guillard et al. 2009). It is un- 
clear how the molecular gas is able to survive the harsh 
ICM environment. Our IRS observations have a reso- 
lution of <~ 3 kpc, insufficient to resolve any structure 
within the molecular hydrogen tail. We do, however, 
see the strongest molecular hydrogen in the part of the 
tail that is closest to the galaxy, and its emission grows 
fainter the further away it is from the galaxy. This likely 
results from both the destruction of the molecular hy- 
drogen by the ICM and the gradual loss of the kinetic 
energy from the ICM-ISM interaction. Given the veloc- 
ity of the galaxy through the cluster, we can place a lower 
limit on the lifetime of the warm molecular hydrogen to 
be 10 — 20 Myr. This number is likely to be larger if 
the molecular hydrogen exists in the Ha and X-ray emit- 
ting parts of the tail, which extend to 40 kpc (Sun et al. 
2007) and 70 kpc (Sun et al. 2006), respectively. How- 
ever, our data do not conclusively support shocks being 
the sole excitation mechanism as X-ray irradiation can 
also be a possible candidate. Further ground-based ob- 
servations characterizing the cold gas (CO and HI mea- 
surements) and hot H 2 gas (near-IR measurements of H 2 
rovibrational lines) are required to fully understand this 
phenomenon. Explaining the existence of warm H 2 in 
the tail from ram-pressure stripping remains a big chal- 
lenge as simulations are only beginning to understand 
the emission properties of the stripped gas tails, which 
do not include any molecular content (e.g. Kapferer et 
al. (2009); Tonnesen & Bryan (2010)). 

The star-forming regions seen outside of the plane of 
the galaxy, some of which are outside the tidal radius, are 
remarkable as a demonstration that ram-pressure strip- 
ping can be a source of intracluster stars. A few of these 
star-forming regions are likely associated with the H 2 
tail indicating that in addition to the warm molecular 
gas there must also be cold molecular gas with sufficient 
density to support star formation. The turbulent flow of 
stripped multiphase ISM of the galaxy may allow for ed- 
dies to form in the wake where star formation can occur. 
This is consistent with ram-pressure simulations that do 
indeed show that significant star formation in high den- 
sity clumps can occur in the wake of the galaxy (Kapferer 
et al. 2009). However, their simulations do not explicitly 
include molecular gas, so the picture painted by the sim- 
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ulation is incomplete and does not fully explain what is 
observed in ESO 137-001. It is entirely reasonable that 
the existence of significant molecular gas in the tail may 
enhance star formation within the tail, but a measure- 
ment of the cold gas fraction is necessary to verify this 
hypothesis. 

Next, we compare the morphology of the tail, the mass 
loss rate, and overall time scale for stripping with simu- 
lations to see if these properties are generally consistent. 
The collimated morphology of ESO 137-001 's tail and its 
moderate mass loss rate indicate that it is in the inter- 
mediate stage of ram-pressure stripping, where the gas 
originally displaced from the gaseous disk of the galaxy is 
permanently removed from the halo of the galaxy (Rocdi- 
ger & Hcnslcr 2005). This phase, which follows the short 
instantaneous stripping phase (<~ 10s Myr), lasts approx- 
imately a few 100 Myr in medium to strong ram-pressure. 
This is broadly consistent with our < 1 Gyr upper limit 
for the ram-pressure stripping timescale. However, these 
simulations do not explicitly include the effects of ram- 
pressure on molecular gas, so it is uncertain if the results 
of simulations can be directly compared. Moreover, it is 
not clear if H2 stripping persists during the entire ram- 
pressure stripping phase. There is already some indica- 
tion that the galaxy may be nearing the end of its ram- 
pressure stripping phase from its lower than expected HI 
gas mass and the small size of the Ha disk within the 
galaxy. 

Our conclusion that the galaxy will lose all of its gas in 
a single pass through the core is also consistent with ex- 
isting observational evidence. Assuming that a warm H 2 
tail is found with either a X-ray or Ha tail, a literature 
survey of galaxies with detected tails in Ha (Gavazzi et 
al. 2001; Cortese et al. 2006; Yagi et al. 2007; Kenney 
et al. 2008; Yoshida et al. 2008; Sun et al. 2007, 2010), 
and/or X-ray (Wang et al. 2004; Sun & Vikhlinin 2005; 
Sun et al. 2006, 2010) reveals that these tails are rare. 
Sun et al. (2006) in their X-ray survey of 25 nearby hot 
clusters only discovered 2 galaxies with X-ray tails, al- 
though another case, ESO 137-002, has now been found 
serendipitously in Abcll 3627 (Sun et al. 2010). How- 
ever, this may not be indicative of the true distribution 
because the X-ray gas in the tail is much cooler than the 
ICM (Sun et al. 2006) and may be impossible to discern 
from the hot ICM X-ray emission especially in the core 
regions where ICM X-ray emission is the highest. More- 
over, the observed Ha tails are low surface-brightness 
features that require relatively long integration times to 
detect them (Gavazzi et al. 2001). We also include the 
case where blue stellar trails were seen in galaxies within 
clusters (Cortese et al. 2007), which is similar to the star- 
forming trail seen in ESO 137-001. Cortese et al. (2007), 
in their HST survey of 13 different clusters at z ~ 0.2, 
only found two spiral galaxies with these types of trails. 
The rarity of these tails implies that the timescale for 
their existence is much shorter than the typical crossing 
time of a massive cluster, which is <~ 5 Gyr. This is 
broadly consistent with our timescale estimate, though 
these features may only be present when stripping is in- 
tense. If warm H2 emission does not directly correspond 
to X-ray or Ha emission or exists for a shorter time, it 
is entirely possible this phenomenon is even more rare. 

5. CONCLUSIONS 



We report on the detection of a warm H2 tail in ESO 
137-001, an infalling cluster spiral in Abell 3627. The 
main results of this study are: 

1 . Through Spitzer IRS spectral mapping observations, 
we detect a warm molecular hydrogen tail (T <~ 130 — 160 
K) extending at least 20 kpc from the galaxy ESO 137- 
001 in the cluster Abell 3627. The full extent of the 
tail is unknown and is limited by our slit coverage. The 
tail is co-aligned with both the X-ray (Sun et al. 2006) 
and the Ha (Sun et al. 2007) tail. We measure a total 
warm Ho mass of 4.1lt°;g x 10 7 M within the H 2 tail. 
Approximately, 50% of the warm H 2 tail mass is located 
outside the tidal radius of the galaxy and is permanently 
stripped. 

2. Multiple temperature components are present in the 
tail, and we fit a two-tempcraturc-component model to 
the observed population distribution. We find that the 
hot component (T = 4721^ K) represents approximately 
1% of the gas mass of the tail in the region that we 
observed. 

3. We identify star-forming regions near the galaxy by 
identifying 8^m excesses (measured relative to 3.6/im), 
which coincide with most of the Ha regions detected by 
Sun et al. (2007). We also identify a large 8^m ex- 
cess centered about the galaxy. Several of these star- 
forming regions exist outside the plane of the galaxy, 
some of which are beyond the tidal radius. Therefore, 
ram-pressure stripping can also be a source of intraclus- 
ter stars. 

4. The molecular hydrogen tail extends farther than 
the extent of the star-forming regions. A comparison 
of the H 2 line luminosity with 24 \im emission strongly 
suggests that star formation is not the main excitation 
mechanism for the warm H 2 tail. 

5. The infrared spectrum of ESO 137-001 itself 
matches expectations for a normal starforming galaxy. 
We estimate a star-forming rate at the nucleus of 0.7 
M yr- 1 . 

6. Assuming a velocity range of 925(cr) < v ga i < 
1602(\/3er) km s _1 for the galaxy, we obtain a warm 
molecular hydrogen mass loss rate of 1.9 — 3.4 M yr -1 . 
We place a strong upper limit on the time it would take 
to strip the remaining gas in the galaxy of 300 — 500 
Myr and one on the total ram-pressure stripping time 
scale of 1 Gyr. This means that the galaxy will have 
lost most of its original gas content after a single pass 
through the cluster core, and that ram-pressure operates 
on timescales much shorter than the crossing time in a 
cluster. These results are consistent with the hypoth- 
esis that ram-pressure stripping can transform galaxies 
quickly when they fall into the cores of massive clusters. 
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Fig. 1. — 6' by 10' Digital Sky Survey Image centered about ESO 137-001 with our SL/LL IRS pointings overlaid. North is up and East 
is left. SL slits are smaller than the LL slits. Regions with both 1st and 2nd order IRS pointings are shown in orange. 1st order only and 
2nd order only pointings are shown in red and blue, respectively. Background subtraction is carried out by outrigger slit pointings. The 
inset is the 2.2' by 4' star-subtracted R^-band image of ESO 137-001 reproduced from Figure 10 of Woudt et al. (2008). This image has 
been contrast enhanced to reveal low surface brightness features. The diagonal line shown in the inset is the direction of elongation of the 
E/S0 population discussed in Woudt et al. (2008), which lies very close in position angle to the X-ray and Ha tails. We chose the long- low 
slits to be coaligned with the X-ray tail discovered by Sun et al. (2006). 
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Fig. 2.— 2' by 2' images of ESO 137-001 at 3.6/wn, Ho (Sun et al. 2007), and rest 17.035^m (H 2 0-0 S(l) transition) wavelengths. 

North is up and East is left. The Ha image is smoothed to 1.2" resolution. In the H2 image, the emission shown is at least 3a above 
the background level. The colors follow the visible spectrum, where blue is represents the faintest emission, and red the brightest. The 
center of the galaxy is shown with a black cross. The polygons outline the spectral extraction regions. The box that is centered about the 
galaxy is termed the nuclear region and includes both short-low (SL) and long-low (LL) coverage. The red polygon includes full SL and 
LL coverage of the tail and is identified as the SL/LL tail region, while the long orange and smaller red rectangular boxes that enclose 
the molecular hydrogen tail have only LL coverage and are termed the LL full tail and far tail regions, respectively. The large green circle 
centered about the galaxy signifies the tidal radius. The 3.6fim image shows a relatively undisturbed spiral galaxy. We observe a cometary 
structure at Slim that strongly suggests ram-pressure stripping. We also see a strong correspondence between the 8^m and Ha images, 
which is expected for star-forming regions. The morphologies at Ha and 8/im do not correspond to the 3.6/xm image, suggesting that the 
star-forming regions are relatively young. A strikingly long molecular hydrogen tail is shown in the 17.035 fim image. The position angles 
of the tails seen in Ha and 8fim are coaligned with the molecular hydrogen tail. This is evidence for warm gas being stripped from a galaxy 
with knots of star formation within it. 
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Fig. 3. — ESO 137-001 LL-only full tail and far tail regions rest-frame spectra. The far tail region flux has been offset by -5 mjy for 
better visibility. The red and blue lines denote the best fit to the full tail and far tail data, respectively. The spectra show strong H2 0-0 
S(l) and S(0) transitions in addition to atomic fine structure lines in both cases. The full tail spectrum reveals some continuum emission 
from warm dust, while continuum emission is virtually non-existent in the far tail case even though strong H2 and fine structure lines are 
still present. This is no surprise because the full tail spectrum contains some of the extraplanar star-forming regions while the far tail 
region was chosen not to include any star-forming regions. 
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Fig. 4. — ESO 137-001 SL/LL tail region rest-frame spectrum. The SL data are shown with a larger stretch within the inset plot. The 
red line denotes the best fit to the data. Atomic fine structure and H2 0-0 S(0) thru S(3), and S(7) lines are detected at greater than 3<j 
significance within this region. The presence of aromatic emission indicates star formation is occurring in this region. 
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Fig. 5. — Excitation diagram for the molecular hydrogen tail. The black, red, and blue points represent the populations in the SL/LL 
tail, LL-only full tail, and LL-only far tail regions, respectively. The LL-only full tail and far tail regions were fit by single temperature 
models, which are represented by the solid red and blue lines, respectively. The SL/LL region population was fit by a two-temperature 
(hot and warm components) model. For the two temperature model fits, the S(0), S(l), S(2), S(3) lines were included. The black curve 
is a two-temperature fit to the SL/LL tail data. The S(7) line is several standard deviations away from the best fit curve. A possible 
solution for this issue may be the inclusion of a third temperature component. However, we have insufficient data to constrain an additional 
temperature component. This is not a significant issue as most of the mass is contained with the lowest temperature component. 




Fig. 6. — 2' by 2' Ha image (left) of ESO 137-001 smoothed to 1.2" resolution and the corresponding stellar continuum subtracted 
Sfim image (right). The colors follow the visible spectrum, where blue is represents the faintest emission, and red the brightest. The 
black cross-hair indicates the center of the galaxy. Sources found by running a source detection algorithm on the Sfim stellar continuum 
subtracted image arc shown as yellow circles in both the Ha and Sfim image. The bright single pixel spots surrounded by a dark annulus 
are residuals from the stellar light subtraction. Every source present in the 8(im image is present in the Ha image with the exception of 
the point source southeast from the galaxy. It is unclear what this source is. The similarity between the two wavelengths establishes the 
Ha and Sfim sources are likely star-forming regions. These sources are in a fan-like configuration downstream of the galaxy. 
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Fig. 7. — ESO 137-001 nuclear region rest-frame spectrum. The red line denotes the best fit to the data. We have full short-low and 
long-low spectrograph coverage in this region. The blue line is a starforming galaxy template spectrum from Smith et al. (2007). The 
template spectrum was rescaled to match the 24 fim flux of ESO 137-001. We see a fairly good match between the two spectra, suggesting 
there is nothing peculiar about the infrared properties of ESO 137-001. 
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Fig. 8.— 2.5' by 2.5' images of ESO 137-001 at 8/xm (left), rest 17.035/im (H 2 0-0 S(l)) wavelength (center), and 24 fan. The colors 
follow the visible spectrum, where blue is represents the faintest emission, and red the brightest. In the 2Afim panel, the orange box is the 
spectral extraction region used to measure the SFR while the blue box is the region used by Sun et al. (2007) to compute their Ha SFR. 
The 8/Ltm data was smoothed to match the I7fim resolution after subtracting the stellar continuum using 3.6/xm data. The center of the 
galaxy is shown with a black cross. While there is a strong correspondence between the 8fim and 2&fim emission, the molecular hydrogen 
tail is much more extended. This hints that the stripped molecular hydrogen must only have sufficient density to form stars shortly after 
it has been removed from the galaxy and that star formation cannot be the main source of excitation. 
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Fig. 9. — A comparison of molecular hydrogen ground state rotational line flux (measured from S(0) thru S(3) transitions) to 24/im flux 
as a function of 24fim flux for SINGS galaxies (Roussel et al. 2007) and warm molecular hydrogen tail. The diamond points represent the 
66 galaxies from the SINGS sample. Two of the significant outliers in the sample, NGC 4450 and NGC 4579, are labelled in the plot. The 
warm H2 in both NGC 4550 and NGC 4579 is thought to be heated through shocks from cloud collisions and/or X-ray irradiation because 
star formation and supernova remnants are unable to explain the high H2 emission. The filled triangular and square points represent 
the values measured in the SL/LL and LL far tail regions of the warm molecular hydrogen tail, respectively. It is clear that the two 
measurements in the tail are significant outliers in this comparison sample and are comparable in value to the ratio for NGC 4450, which 
is likely heated by shocks through cloud collisions in addition to X-ray irradiation (Roussel et al. 2007). 
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Fig. 10. — Ratio of cluster radial a r and galaxy centripetal a ga i acceleration (a r /a ga i) and effectiveness of ram-pressure ((jyam) as a 
function of cluster-centric distance. The solid and short dashed curves represent a r /a ga l determined from a best-fit NFW cluster mass 
profile and a X-ray-data-dcrivcd mass profile, respectively. If the ratio is greater than 1 (the dotted horizontal line), the cluster tidal force 
will be sufficient to tidally strip the galaxy. However, at the projected distance of 280 kpc (the dotted vertical line), the minimum distance 
the cluster galaxy can be from the center, the a r /a ga i is less than 1 for both cluster mass profiles, indicating that tidal stripping from 
the cluster gravitational field is not a significant effect. The dash dotted and long dashed curves represent q rarn for a galaxy traveling at 
925 and 1602 km s _1 , respectively. When q ra m > 1, the Gunn & Gott (1972) criterion is satisfied and ram-pressure stripping becomes 
effective. It is clear that ram-pressure stripping dominates over tidal effects and begins at fairly large cluster centric distances: 720 kpc 
(925 km s- 1 ) and 1.4 Mpc (1602 km s" 1 ). 
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TABLE 1 

Measured H2 Rotational Line Fluxes for ESO 137-001 Tail 



Spectral Feature 


A 


Tail (SL/LL) Flux 


All Tail (LL-only) Flux 


Far Tail (LL-only) Flux 




W 


(10- 17 W m" 2 ) 


(10- 17 W m- 2 ) 


(IO" 17 W m" 2 ) 


0-0 S(7) 


5.511 


1.03±0.26 






0-0 S(6) 


6.109 


< 0.74* 






0-0 S(5) 


6.909 


< 1.04* 






0-0 S(4) 


8.026 


< 0.59* 






0-0 S(3) 


9.665 


1.20 ±0.11 






0-0 S(2) 


12.279 


0.42 ±0.08 






0-0 S(l) 


17.035 


0.92 ±0.01 


2.10 ±0.08 


0.80 ±0.09 


0-0 S(0) 


28.221 


0.23 ±0.01 


0.47 ±0.03 


0.21 ±0.02 



3<t detection limit. 



TABLE 2 

Fine Structure Line Fluxes for ESO 137-001 Tail 



Spectral Feature A Tail (SL/LL) Flux Full Tail (LL-only) Flux Far Tail (LL-only) Flux 

{urn) (10- 17 W m" 2 ) (IO" 17 W m" 2 ) (IO" 17 W m" 2 ) 



Nell] 12.8 1.17 ±0.11 

Nelll] 15.6 0.25 ±0.01 0.17 ±0.08 < 0.21 

SIII] 18.7 0.30 ±0.01 0.34 ±0.08 0.17 ±0.04 

OIV],[FeII]t 25.9,26.0 0.14 ±0.01 < 0.10* 

SIII] 33.5 0.87 ±0.02 0.90 ±0.10 0.18 ±0.06 

Sill] 3JL8 1.16 ± 0.02 1.35 ±0.11 0.57 ± 0.08 



t Possible blended feature. 
3<r detection limit. 



TABLE 3 

Measured H 2 Gas Masses for ESO 137-001 



Fit 


Component 


Tex (K) 


Mot (cm" 2 ) 


S (M© pc- 2 ) 


H 2 Mass (M ) 


Full Tail (LL-only) 


Warm 


157 ±3 


2.56t°; 2 J x 1019 


41+ - 04 


1.1lt° f 9 x 10 7 


Far Tail (LL-only) 


Warm 


149 ±6 


2.9O+0 43 x 10 19 


0.46_ 007 


2 - 15 -0%2 X 107 


Tail (SL/LL) 


Warm 
Hot 


470+ 77 


5.9+j- 2 x 10 19 
5.9l^ x 10 17 


Q4+ u - bY 
9At 7 4 ° x IO" 3 


3.6+0^ x IO 7 
3.6l 2 'g x 10 5 



TABLE 4 

Dusty 8/j.m Source List 



# 


RA 


DEC 


Distance (") 


Ha counterpart 


1 


1G.1 5.07 in 

16:13:2 { .19 


-uU:4o:oU.44 


U.8 


Y 


2 


i/" 1 io O/i no 

16:13:24.98 


-oU:45:43.51 


18.3 


Y 


3 


16:13:26.48 


(~T\ /I ET OO OCT 

-00:45:32.85 


18.6 


Y 


1 


16:13:25.45 


-0U:45:o4. /9 


20.7 


Y 


5 


16:13:23.93 


-60:45:53.60 


24.9 


* 

Y 


6 


16:13:24.54 


-60:46:05.92 


25.4 


Y 


7 


16:13:24.13 


-60:45:33.46 


28.8 


Y 


8 


16:13:23.66 


-60:45:23.70 


37.8 


Y 


9 


16:13:25.78 


-60:45:10.21 


41.8 


Y 


10 


16:13:29.82 


-60:46:30.42 


44.0 


N 


11 


16:13:22.83 


-60:45:20.55 


44.4 


Y 


12 


16:13:23.03 


-60:44:49.45 


68.6 




13 


16:13:35.89 


-60:45:23.07 


68.6 


Y 


14 


16:13:33.74 


-60:44:49.60 


77.0 


N 


15 


16:13:32.01 


-60:47:00.30 


77.9 


N 


16 


16:13:36.16 


-60:44:53.67 


86.2 




17 


16:13:32.05 


-60:44:27.37 


90.1 


N 


18 


16:13:30.07 


-60:47:25.82 


97.5 




19 


16:13:13.65 


-60:46:04.68 


101.0 


N 


20 


16:13:39.49 


-60:44:47.56 


109.3 


N 


21 


16:13:42.58 


-60:46:28.59 


118.3 





This single source detection may actually be two separate 
sources. 



TABLE 5 
Aromatic Feature Fluxes 



Aromatic Feature 


Nuclear Flux 


Tail (SL/LL Region) Flux 




(10- 17 W m" 2 ) 


(10~ 17 W m" 2 ) 


6.2 (im 


66.1 ± 1.8 


2.21 ±0.72 


7.7 /jm Complex 


217 ± 7 


21.0 ± 2.1 


11.3 fim Complex 


57.8 ±0.9 


3.43 ± 0.29 


12.6 fim Complex 


32.9 ±0.9 


0.88 ±0.35 


17.7 (im Complex 


28.1 ±0.3 


1.76 ±0.07 



